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Abstract 

We generalize the non-linear one-dimensional equation of a fluid layer for 
any depth and length as an infinite order diff'erential equation for the steady 
waves. This equation can be written as a q-differential one, with its general 
solution written as a power series expansion with coefficients satisfying a nonlinear 
recurrence relation. In the limit of long and shallow water (shallow channels) we 
reobtain the well known KdV equation together with its single-soliton solution. 
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1 Introduction 



The Korteweg-de Vries equation (KdV) was proposed by Korteweg and de Vries [1] 
one century ago to explain the steady translation water waves observed in a channel 
by Scott Russel in 1834. The main property of this equation consists in the equal 
occurence of non-linearity and dispersion. Some of its most important solutions, the 
solitary waves, are single, localised and non-dispersive entities (no dispersion in shape), 
having also locahsed finite energy density. Among these solutions, the solitons are 
solitary waves with special added requirements concerning their behaviour at infinity 
(x — > oo, t ^ oo) and having special properties concerning the scattering between 
different such single-solitary solutions. Despite this non-hnearity, the KdV equation 
is an infinite dimensional Hamiltonian system [2-5] and, remarkable, when the KdV 
solution evolves in time, the eigenvalues of the associated Sturm-Liouville operator 
d'^ + r], where rj is the solution of the KdV equation, remain constant. Several extended 
and complex methods have been developed to study and to solve the KdV equation and 
other non-linear wave equations (non-linear Schrodinger equation, NLS, modified-KdV 
equation, MKdV, sine-Gordon equation, etc.): inverse scattering method, theoretical 
group methods, etc. (for a review of some of these techniques one can look into the 
book by Witham [6] or the paper by Scott ct al [7] . 

Due to its properties, the KdV equation was the source of many applications and 
results in a large area of non-linear physics (For a recent review see [8] and the references 
herein). For example, the KdV equation and its generalizations (KdV hierarchies, KP 
equation and its hierachies and their supersymmetric generalizations [9]) appear in 
different fields of physics, from plasma physics [10] to string theory [11], or to describe 
the dynamics of quantum extended objects [12], and nonlinear phenomena in nuclear 
physics (non-linear liquid drop model) [13,14], and the list is far from being exhausted. 
In all these applications, the KdV equation (or other non- linear ones) appears as the 
consequence of certain simplifications of the physical models, especially due to some 
perturbation techniques involved. 
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In the present paper we have used only one of the two well-known necessary 
conditions for obtaining the KdV equation in shallow water channels, i.e. the smallness 
of the amphtude of the sohton with respect to the depth of the channel, h. This is the 
single condition which we use, and hence, we consider h to be an arbitrar parameter. 
In order to obtain the KdV equation from shallow water channel models [3,5], a second 
condition is imposed, i.e. the depth h of the channel should be smaller than the half- 
width of the solitary wave. We do not use, in our generalised approach, nor this second 
condition neither the condition of infinite length of the channel. In our case the length 
of the channel, L, is also considered an arbitrary parameter. 

We introduce this generalisation in connection with the last (above mentioned) nu- 
clear physics application. If the height of the perturbation is comparable to the nuclear 
radius (big clusters formation on the nuclear surface or symmetric nuclear molecules) 
it is not appropriate to consider that the nuclear dynamics is that of a shallow fluid 
layer. Therefore, we should like to study, as a flrst step of our aim, the non-linear dy- 
namics of a fluid in a channel of arbitrary depth, in a bounded domain. This different 
starting point leads us to a new type of equation which generalizes in some sense the 
KdV equation and reduces to this in the shallow water case. 

As we shall show, this generalized equation (inflnite order differential equation) 
can also be written in terms of a q-differential equation. In this context we introduce 
here another possible physical interpretation for the q-deformation parameter, in the 
fleld of fluid dynamics, as the depth of the fluid layer. We use the connection between 
the theory of quantum algebras (or q-groups, q-deformed algebras), originally initiated 
by Jimbo and Drinfeld [15], and the q-differential equations [16]. 

In the present model we take into account not only the gravity but also the 
influence of the surface pressure acting on the free surface of the fluid. This imphes 
that the KdV-like structure for the dynamical equation arrises in a lower order (first 
order in the smallness parameters) and one may not take into account the second order 
of smallness to obtain a KdV equation. We also note that, despite the general tendency 
followed by papers concerning applications of q-groups in physics of flrst constructing 
the deformed algebra and later searching possible applications, we have obtained a 
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q-differential equation in a natural way, by starting from a traditional physical one- 
dimensional hydrodynamical model. 



2 The generalized KdV equation 



Let us consider a one-dimensional ideal incompressible fluid layer with depth h and 
constant density p, in an uniform gravitational field. We suppose irrotational motion 
and consequently the field of velocities is obtained from a potential function y, t) 

V{x,y,t) = {u,v) = V^x,y,t). (1) 

The continuity equation for the fiuid results in the Laplace equation for 

A^x,y,t) = 0. (2) 

The above Laplace equation should be solved with appropriate boundary conditions 
for our physical problem. We take a two-dimensional domain: x e [xq — L, xq + L] 
(the "horizontal" coordinate) and y G [0,.^(a;, t)] (the "vertical" coordinate) where Xq 
is an arbitrar parameter, L is an arbitrary lenght and ^{x, t) is the equation of the free 
surface of the fluid. The boundary conditions on the lateral walls, x — xo±L and on the 
bottom y — 0, consist in the condition of vanishing of the normal velocity component. 
The free surface fulfills the kinematic condition [3] : 



(3) 



s 

where we denote by E the free surface of equation y — ^{x,t) and the lower index 
indicates the derivative. Eq. (3) expresses the fact that the fluid particles, which 
belong to the surface, remain durring time evolution in this surface. By taking into 
account the above boundary conditions on the lateral walls and on the bottom, we can 
write the potential of the velocities in the form: 

^x,y,t) = f;cosh(^)(a,(t)cos(^^) +/3fc(t)sin(^^)), (4) 



where ak and (3k are time dependent coefficients fulfiling tfie condition 

. kn(xo ^ L) kirixQ ± L) 
ctfesm = Pfecos , 

for any positive integer k. Tfiis restriction introduces a special time dependence of 
and i.e. ^ = 7fc = const. ^ for any A;, when f3k 7^ 0. If = we simply equate with 
the inverse of the above fraction. It exists also a dependence of the ak,f3k functions 
of Xq. This special time dependence does not affect the generality of the potential 
but only the balance between the two terms in the RHS of eq.(4). If we fix all Pk — ^ 
we can take arbitrary values for xq and general form for ak- In the infinite channel 
limit, L — > 00, there are no more restrictions concerning and (3k functions, and we 
can choose xq = 0, too. 

We introduce the function: 

fi^^t) = ^ ^(^-Qfe(t)sin(^A;7r^^ +/3fc(t)cos(^/c7r^^^, (5) 

and the velocity field can be written like: 

u = ^ cos(yd)f(x,t) 

V = %^-sm{yd)f{x,t). (6) 

where, for simplicity, the operator d represents, from now on, the partial derivative with 
respect to the x coordinate. The equations (5,6) do not depend on L and we equally 
treat the case L — > 00, when we look for unbounded domains and travelling profiles for 

the solutions. 

In the following we use in the boundary condition (3) the velocities evaluated at 
y = C,{x,t) = h + r){x,t) in the first order in 77, as we describe the perturbations with 
small height compared to the depth, but not necessarly with large widths. Eqs.(6) read: 

u{x,^{x,t),t) = [cos{hd) — ri{x,t)dsm{hd)] f{x,t) 

v{x,^{x,t),t) = -[sm{hd)+r]{x,t)dcos{hd)]f{x,t). (7) 

The dynamics of the fluid is described by the Euler equation at the free surface. 
The Euler equation, written on E, results (after writting it in terms of the potential, 
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differentiating it witli respect to x and imposing tfie condition y — $^{x,t)) in tlie form 

1 

ut + uux + vvx + gVx + -Px = 0, (8) 

wfiere g represents the gravitational acceleration and P the surface pressure. Following 
the same approach used in the calculation of surface capilar waves [17], we have for our 
one-dimensional case 



P 



= ^ = — '^^^^ , /o — —cVxx, for small ri. (9) 



where TZ is the local radius of curvature of the surface (in this case, of the curve 
y — ^{x,t)) and a is the surface pressure coefficient. In the interior of the fluid the 
pressure is given by the Euler equation. The nonlinearities appear in the dynamics 
through the nonlinear terms in eqs. (3), (7), (8) and (9). Consequently, we have a 
system of two differential equations (3,8) for the two unknown functions: f{x,t) and 
r]{x, t), with u and v depending on r] and / from eqs. (7). With / and r] determined and 
introduced in the expressions of u and v, we can finally find the coefficients and Pk- 
In the following we treat these equations in the approximation of small perturbations 
of the surface E, with respect to the depth, a — max\r]^'^\x, t)\ « h, where k — Q, 3 
are orders of differentiation. 

In the linear approximation eqs. (3) and (8) become, respectively 



-sm{hd)f = r)t 



a 



cos{hd)ft = -gr]x + -Vxxx: (10) 

P 

and we obtain, by eliminating from eqs. (10), 

cos{hd)r)tt = sm{hd)(J^r)^ (11) 

which, in the lowest order of approximation in (hd) for the sine and cosine functions, and 

in the absence of the surface pressure, gives us the familiar wave equation rju = c^rixx, 
where Cq = ^/gh is the sound wave velocity. By introducing the solution r] = e'^i^^-^^) 
in the linearized eq.(ll) we obtain a non-linear dispersion relation 

a \ tanh(A;/i) 



a;2 



V pgJ 



pgj kh 
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In the limit of shallow waters we find, for the dispersion relations, in the case a — 
the acoustic waves limit, and, for a >> pg the surface capilar waves limit uj"^ — ^k^- 
In the absence of the surface pressure {a — 0) the function / is given, in this linear 
approximation, at least formally, by: 



= ^(^^^^^7^) V, (12) 



Co f sm{2hd) y'^^ 
J V 2hd J 



which in the limit of shallow fluid has a particular solution in the form 

For the time derivative of / we have, from the second equation of eqs.(lO),: 

d cr 



fl'^ix,t) = (cos{hd))-'{-^Vx + -Vx..), (13) 



which in the limit of a shallow fluid reduces to: 



n p 



XXX 



Following [3], we look for the solution of eqs.(3,8) in the form 



/ = j^^+y/. ^ (14) 

ft = -Co^(cos(/ia))"^77a; + ^(cos(/ia))"^?7xxa; + (^^^ 5f2, (15) 

which represents a sort of perturbation technique in a/h, where rj = afj. Of course a 
functional connection exists between the perturbations f2{x,t) and g2{x,t). Eq.(3), in 
the lowest order in a/h, yields: 

— Cosm{hd)fj = hfjt + acQ{fjxCos{hd)fj + fjcos{hd)fix). (16) 

If we approximate sin(/i9) ~ /i9 — |(/i9)^, cos{hd) ~ 1 — |(/i9)^, we obtain, from 
eq.(16), the polynomial differential equation: 

^ 2, ~~ ,~ ~ coe^h^ f ^ _ \ r. 

ar]t + 2coe hr]ri:, + coehr]:c - coe—rj^^^ — hlxVxx + VVxxxj = 0, (17) 

where e = f • The first four terms in eq.(17) correspond to the zero order approximation 
terms in /, obtained from the boundary condition at the free surface, in eq.(6.1.15 a) 



from [3], i.e. the traditional way of obtaining the KdV equation in shallow channels. 
In this case all the terms are in first and second order in e. If we apply the second 
restriction with respect to the solutions, i.e. the half-width to be much larger than h, 
we can neglect in eq.(17) the last paranthesis, and we obtain exactly the KdV equation 
for the free surface boundary condition. In other words we understand the condition 
hd "small" like {hd)f{x,t) <^ 1 all over the domain of definition of /. This means that 
the spatial extension of the perturbation f{x, t) is large compared to h, which is exactly 
the case in which KdV equation arrises from the shallow water model (see Chapter 6 
in [3], hdf(x, t) of order ~ | = 5 < 1). 

By using again the approximations given by eq.(7), we can write the Euler eq.(8) 
in the form: 

dtnf + nf{dtnf) + nf{dnf)+u;f{du;f) = -gv^ + ^^^x, (18) 
where we used the notations: 

fl — cos{hd) — rid sm{hd) 

iij — sm{hd) + r]dcos{hd). (19) 
We note here an interesting property of the operators given in eq.(19): 



fl + iu} = e 



iyd 



02{hd) + 02{hri). (20) 



s 

Following the same procedure like for the free surface boundary condition eq.(3), i.e. 
using the approximation for small r^, we obtain, from eq.(18): 

cos{hd)ft = -grjx + -Vxxx, (21) 

P 

which, in the lowest order in a/h and by using eqs. (14,15), reduces to an identity. 

Before further analysis, we would like to note that in the shallow water case, 
following the same addimensional notations like in Chapter 6 of [3], i.e. e = ^, S — j, 
where I gives the order of magnitude of the half-width of the perturbation 77, and 
introducing the new addimensional parameter a — we obtain, for the Euler 

equation, the form: 

Vt' + Vx' + ^VVx' + a^Vx'x' = 0, (22) 



which is again the traditional KdV equation. The primes atached to the lower indexes 
represent addimensional units, according to [3]. The difference between eq.(22) and the 
corresponding eq.(6.1.15.b) from [3] is given by the inclusion of the surface pressure 
effects. If the coefficient a. is large enough, one is no more obliged to take into account 

second order terms, i.e. 5^, to obtain the KdV equation. Of course, the above approach 
introduces changes in the differential equations which are involved with higher order 
perturbations, i.e. /^^^ and /*^^^ from [3] or /2, gi in our case. 

The reduction of eq.(3) to the KdV equation occures if, in eq.(16) we limit only 
to third order of magnitude in hd: 

^2co „ 
6 h 

In the following we shall investigate this generalised KdV equation (GKdV), obtained 
from eq.(16), by keeping all terms in sin and cos, namely: 

+ ^ sin(/i9)?7 + ^{j]x cos{hd)7] + 7] cos{hd)7]x) = (23) 

ft ft 

Eqs.(3) and (18), in higher orders in | and in hd, yield the corresponding differential 
equations for the functions /2 and g2, but we shall study here only eq.(23). 



3 The q-differential form of the GKdV equation 



In this Section we limit ourselves to the steady translation waves and consider only 

solutions of the form i]{x,t) = rj^x + Acot) = f]{X) where A & R and X = x + Acot. 
Eq.(23) can be written in the form: 

2i 2 
^^^^^ryx(^+|M±Z^x(X-^ = 0, (24) 
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if we suppose that 77 is an analj^ic function in the domain Re{z) e (00,00), Im{z) e 
{—h,h) of the complex plane. We study the solutions with a rapid decreas at infinity 
and make a change of variable: v — e^^ for x e (—00, 0) and v — e~^^ for x e (0, 00). 
Here S is a positive real constant. Putting: 

rj{X) = -hA + f{v), 

we obtain a q-differential equation for the function f{v): 



(25) 



DMv) ' DJ(v) ■ - B 

where we denoted q = e*^'*, {\q\ = 1), and the operator of the q-derivative is defined 
[16]: 

f{qr)-fiq-'r) 



qv — q '^v 

We then write the solution of eq.(24) (or cq.(25)) as a power series in v: 

00 

f(v) = $^a„t;", 

n=0 



(26) 



(27) 



and we must put oq = hA to have limx-^±car]{x) — 0. Equation (25) results in a 
recurrence non-linear relation for the coefficients a„: 



Ahk + — sm{Bhk) 
B 



k-i , . 
= ~ 51 ^( {Bh{k — n)) + cos{Bh{k — 1)) ja„ajfc_„. 
n=l ^ ^ 



(28) 



By taking k — 1 in the above relation we obtain oi Ah + ^sin(S/i)j = 0. Without 
any loss of generality, due to the arbitraryness of B we can write: 

sm{Bh) 



A = 



Bh 



(29) 



This relation fixes the velocity of the envelope of the perturbation if its asymptotic 
behaviour is fulfilled. To have A 7^ we need Bh ^ kn for k integer. In this condition 
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oi is arbitrary at present and, by writing = ccfeOi we have ai — 1 and the recurrence 
relation: 

25 cos ^^^^ fe-i Bh(2k-n-l) , , 

= A.sin(5/.)-sin(fc5/.)£"^°^ 2 ^'^^ 

for k > 2. This recurrence relation gives the A; coefficient if we know the first k — 1 
ones. 

To obtain a smooth behaviour of the solution r]{X) in X = 0, i.e. a continuity of 
its derivative, we must introduce the condition: 

oo 

/,(!) = ^naX"' = 0, (31) 

n=l 

or the derivative of the power series f{v) with the coefficients given by eq.(28) is zero 
in z & R, z — ai. This condition fixes oi. 

We study now a limiting case of the relation (28), by replacing sin and cos ex- 
pressions with their lowest nonvanishing terms in their power expansions. Thus, we 
have: 

6 

^ i^2;,3^(;t2_i)E^«n«^- (32) 
It is a straightforward exercise to prove that: 



1 



k-l 



is the solution of the above recurrence relation. This can be done by using the mathe- 
matical induction and by taking into account the relations: 

2 _ A;(A;-1)(2A;-1) 



E 



n — 



We can write the power expansion: 

oo 



oo ^ 1 s k-l 

9i^) = EM^T^ (34) 



k=l 
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having the radius of convergence TZ — 2B'^h^ (due to the Cauchy-Hadamard criterium). 
The function g(z) can be written in the form 



1 



2B2/i3 

The condition (31) results in oi = —2B'^h^ and 



^ ^ 2B^h' = 2 (35) 

(1- 



= H^y~\-^B'h'r = 2B'h'{-l)% (36) 



which gives: 



r){X) = 2^2/^3 g A; (-e-^l^l 



p-m\ 

2B^h^ 



a + e-^W)' 
E^fe3 1 

2 (cosh(5X/2))^' 



(37) 



As one expects, the above solution is exactly the single-soliton solution of the KdV 
equation and it was indeed obtained by assuming h small in the recurrence relation (28). 
We stress at this point that the GKdV equation has two general features, expressed 
in the reduction of eq.(16) to KdV equation and in eq.(23): in the hmit hd "small", 
both the differential equation GKdV and one of its solutions, r]{X), go versus the KdV 
equation and its single - soliton solution, correspondingly. 

In the general case we do not obtain a simple solution like the relation (33) but 
we have all the necessary informations from the relations (28). It seems that the power 
series g{z) with the coefficients given by the recurrence relations (28) has a nonvanishing 
radius of convergence. The problem of the existence of a real point Zq in the disk of 
convergence or on its frontier in which ^''(^o) = needs further studies. For the KdV 
equation this point is on the frontier of the disk of convergence of the power series with 
coefficients given by the relations (28). We mention that g'{zo) — means that the 
function is not univalent ( it is not injective in the vicinity of this point). 

An interesting comment concernig an unexpected connection is appropriate: Bieber- 
bach formulated in 1916 the following conjecture: if f{z) = z + a2z'^ + ... is analytic 
and univalent in the unit disk, then |a„| < n for all n, with equality occuring only for 
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rotations of the Koebe function k{z) — jj^^ [18]. The form of our solution for the KdV 
equation, eq. (34), is exactly the Koebe function up to a scaling of the z variable. This 
conjecture was demonstrated in 1986 by de Branges [19]. This fact suggests that one 
could obtain information about the possible zeros of the g'{z) if it is possible to obtain 
estimates for the convergence radius of the power series with the coefficients given by 
(30). 

The change of the variable which we used to obtain the non-linear q-differential 
equation (25) can be used in the time-dependent equation obtained from eq.(23). If 
tjIx, t) — f{v, t) where v — e^^ or v — e~^^, then we obtain: 

CqV Dq}^{V,t) 

±Briv)^^^^j^±2sm{Bh)D,r{v,t) = 0, (38) 

where ± reffers to the positive or negative real semiaxis. If we try a solution of the 
form f^{v,t) — Z^^o^fe'"'^) then we obtain the following infinite system of differential 
equations: 

Co at 

Bh(n-l)J^ , ^ /(2k-n-l)Bh\ , , 

zp 2i3cos^^ lY.^nal{t)al_^{t)^o^\^- ^-^J- (39) 

Wc note that if ao(^o) = then this relation is preserved in time. The above system 
can be resolved step by step as the equation for a„ involves only the coefficients a/j with 
k < n. However, the smooth behaviour of the solution in a; = at every moment is a 
nontrivial problem. 



4 Conclusions 



We have shown that the KdV equation describing the shallow liquids can be generalised 
for any depths and lengths. Consequently the influence of the layer depth could be 
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included in the present studies on different fields, like clusterization on the nuclear 
surface in nuclear physics. 

The present paper shows two important results: a generalization of the KdV 
equation starting from a physical model and an embeding of this non-linear equation 
into a q-dcformcd formalism. 

By using a non-linear hydrodynamic approach for a fluid layer wc obtained a 
dynamical differential equation, of infinite order, wich generalizes the KdV equation 
for shallow water. Contrary to the later one, the generalised equation is valid for any 
depth and length of the fluid layer. We succeded to write this equation in a q-differential 
formalism, and consequently we have obtained a non-linear recurrence relation for the 
coefficients of its general solution. We stress the importance of the introduction of the 
surface pressure term, eqs.(8,9), which provides the necessary dispersion term, i.e. r]xxx- 
This term is esential in two respects: first, it introduces the dispersion in a smaller order 
than in the traditional case [3], and second, it represents the single term responsible 
for dispersion in the case of cyhndrical geometry [13,14]. Both the generahzed KdV 
equation and its formal solution reach the KdV system in the shallow water limiting 
case, which gives confidence in the present approach. We conjecture that exists a deep 
connection between non-linear differential equations, infinite order linear differential 
equations (which are connected with finite diference equations) and q-diffcrcntial equa- 
tions and their q-dcformed symmetries. In this sens, in the present paper we succeded 
to show that by starting from a one dimensional model for an ideal incompressible ir- 
rotational fluid layer, a more general differential equation, than the KdV equation, was 
obtained. The equation which we proposed not only generalizes the KdV equation, but 
also can be written as a q-deformed nonlinear equation with a deformation parameter 
of modulus one. In our case the deformation parameter q is put into relation with the 
depth of the fluid. To our knowledge this one of the few cases when one can add a 
physical interpretation to this parameter. 

We stress that the result of the Section 3, i.e. the realisation of the GKdV equa- 
tion as a q-deformed one, could be the starting point to search for more interesting 
symmetries. Perhaps, there is a more general connection with the q-deformed algebras. 
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It would be interesting to interpret the q-GKdV eqs. (24,25) as the Casimir element of a 
certain q-deformed algebra. This posibility could open the way to construct invariants 
of this equation and to compare them, in the limiting case of shallow waters, with the 
KdV invariants. Such work is in progress. 

We expect many possible application of such a formalism, especially in the field 
of nuclear clusterization. If the GKdV equation arrises from a Lagrangean formalism 
or from a Hamiltonian formalism, then, one can apply this result in physical models 
which involve exotic shapes (cluster decay, spontaneous fission, fragmentation of nuclei, 
atomic clusters, etc). 
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